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FIELD OF THE INVENTION 

This invention relates to a method for classifying heart rhythms based on the shape (morphology) of the 
endocardial or epicardial electrogram. More particularly, a number of features (metrics) of the electrogram are 
5 measured to form a description of the shape of individual electrograms, and the metrics are then converted 
into a cardiac rhythm diagnosis by means of a pattern classification technique. 

BACKGROUND OF THE INVENTION 

10 There are a number of ways to discriminate different types of heart rhythm in an implantable cardioverter 
defibrillator (ICD). The most obvious way to do this is to use heart rate information alone, in the manner shown 
for example in U.S. Patent No. 4,475,551 to Langeret al., entitled "Arrhythmia Detection and Defibrillation Sys- 
tem and Method", which issued October 9, 1984. The purpose of discriminating different rhythms in an ICD 
device is to allow the device to respond appropriately to each of them. Problems occur if there are two or more 

15 rhythms which require a different response from the device but cannot easily be discriminated from one an- 
other. 

The rhythms which cause most problems in existing ICD devices are atrial fibrillation (AF) and sinus ta- 
chycardia (ST). Both of these rhythms can result in a ventricular rate that is high enough for it to be considered 
to be ventricular tachycardia (VT). Thus, VT cannot be discriminated from AF or ST on the basis of ventricular 

20 rate alone. Additional information must be used to correctly identify these rhythms. 

As discussed in an article by Arzbaecher et al., entitled "Automatic Tachycardia Recognition", appearing 
in PACE . May-June 1 984, Volume 7 (II), pages 541-547, and in an article by Jenkins et al., entitled "Tachycardia 
Detection in Implantable Antitachycardia Devices", appearing in PACE , Volume 7 (II), pages 1273-1277, No- 
vember-December 1984, atria! fibrillation can be readily identified by looking at the ratio of the atrial rate to 

25 the ventricular rate. This requires the use of an atria! lead, however, which may be undesirable or unavailable. 
Also, this method is unsuited to discrimination of ST and VT, since both rhythms may exhibit a 1 :1 ratio of atrial 
to ventricular rate. 

Timing methods have been proposed for distinguishing ST or AF from VT (see, e.g., U.S. Patent No. 
4,830,006 to Haluska eta!., entitled "Implantable Cardiac Stimulatorfor Detection and Treatment of Ventricular 
30 Arrhythmias", which issued on May 16, 1989, and see, e.g., an article by Camm et al., entitled Tachycardia 
Recognition by Implantable Electronic Devices", appearing in PACE , September-October 1987, volume 10, pa- 
ges 1175-1190). These methods are of dubious accuracy and some of them also require the use of an atrial 
lead. 

The shape or morphology of the intracardiac electrogram can be used to achieve the desired discrimination 

35 between ST and VT. 

Template methods of morphology analysis provide a scalar measure of the difference between two shapes. 
An example of normal morphology is established, and this is compared to each heartbeat for analysis. Ascalar 
bound or threshold on the normal variation in this scalar measure is established, and any measure that exceeds 
this threshold is considered to represent an abnormal beat. 

40 A number of template methods have been heretofore described (see, e.g., U.S. Patent No. 5,000,189 to 
Throne et al., entitled "Method and System for Monitoring Electrocardiographic Signals and Detecting a Path- 
ological Cardiac Arrhythmia Such as Ventricular Tachycardia", which issued on March 19, 1991, and see, e.g., 
an article by Throne et al., entitled "A Comparison of Four New Time-Domain Techniques for Discriminating 
Monomorphic Ventricular Tachycardia from Sinus Rhythm Using Ventricular Waveform Morphology", appear- 

45 'ng in IEEE Tr. on BioMed. Eng .. Vol. 38 (6), at pages 561-570, (June 1991)). Template methods in general 
require high sampling rates and high processing overheads. For these reasons they are unsuited for use in an 
ICD. 

Many attempts have been made to find a single metric that can be used for rhythm discrimination purposes 
(see, e.g., an article by Pannizzo et al., entitled "Automatic Methods for Detection of Tachyarrhythmias by An- 

50 titachycardia Devices", appearing in the Journal of American College of Cardiology at Vol. 1 1 , pages 308-31 6, 
February 1 988). A probability density function (PDF) algorithm has been of some utility in diagnosing ventric- 
ular fibrillation (VF), but has not proven successful for discriminating ST from VT. The PDF algorithm is dis- 
cussed in the aforementioned U.S. Patent No. 4,475,551 to Langer et a!., and in an article by Bach et al., entitled 
"Engineering Aspects of Implantable Defibrillators," appearing in a book edited by Saksena et al., entitled "Elec- 

55 trical Therapy for Cardiac Arrhythmias: Pacing, Antitachycardia Devices, Catheter Ablation", published by W. 
B. Saunders, Philadelphia, 1 990, at pages 375-376. 

U.S. Patent No. 5,086,772 to Larnard et al. describes a method for combining two simple morphological 
features with timing information to improve the rhythm discrimination process. This method considers first the 
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rate, and uses the morphological information only in a specified rate band, to classify individual heartbeats. 
The cardiac rhythm is then diagnosed on the basis of the classification of a number of successive heartbeats. 
The present invention is a generalization of the Larnard et al. morphological method but its application is not 
limited thereby. 

5 It is, therefore, a primary object of this invention to provide an improved method of classifying heart 

rhythms by means of electrogram morphology. 

It is another object of this invention to provide an improved method of classifying heart rhythms by means 
of electrogram morphology, which method is usable in an implantable cardioverter defibrillator. 

It is a still further object of this invention to provide a method of classifying heart rhythms by means of 
10 electrogram morphology, which method employs a low sampling rate, requires very few multiplications per 
heartbeat classification, consumes very little power in the classification process, and can be implemented 
both in digital and analog electronics. 

Further objects, features and advantages of the invention will become apparent as the following descrip- 
tion proceeds. 

15 

SUMMARY OF THE INVENTION 

Briefly stated, and in accordance with one embodiment of this invention, a number of metrics are used 
simultaneously to improve the accuracy of the classification decision process. Bounds on the normal variation 

20 of the metrics are established. These bounds define a decision surface which divides the measurement space 
into two regions. Measurement vectors which fall in one of these regions are considered to represent normal 
heartbeats, while those that fall in the other region are considered to represent abnormal heartbeats. Methods 
for achieving this have been discussed by K. Fukunaga in a book entitled "Introduction to Statistical Pattern 
Recognition", published by Academic Press, Inc., London, Second Edition, pp. 124-169, 1990, and by Duda 

25 et al. f in a book entitled "Pattern Classification and Scene Analysis", published by John Wiley & Sons, New 
York, pp. 24-31, 1973. The particular metrics disclosed in the present invention require only simple operations 
and low sample rates to calculate. Further, the method for reaching a rhythm classification based on these 
metrics imposes a very small computational burden. Thus, this method is ideally suited for use in ICDs. 
The shape of the electrogram is quantified by measuring particular features thereof, such as amplitude 

30 or width. The value of these features are the metrics, and a set of metrics form a quantitative description of 
the shape of an electrogram. The shape of an electrogram is then indirectly compared to that of a known normal 
electrogram by comparing the metrics, and a given heartbeat is then classified as either normal or abnormal. 
VT may then be diagnosed if a number of successive heartbeats are abnormal. 

The metrics that are discussed here are simple ones that can be implemented either in analog electronics, 

35 or with low computional overhead on a microprocessor. They require only such techniques as filtering, peak- 
picking, thresholding, and integrating. Further, by use of a peak-holding sampling process, they can be calcu- 
lated from a very low data rate digital signal. The rate of 80 Hz has been shown to be sufficient A suitable set 
of metrics that may be employed in connection with diagnosing a rhythm are the area under the electrogram 
(hereinafter defined more completely and referred to as "MOD"), the peak-to-peak amplitude (hereinafter de- 

40 fined more completely and referred to as "PKP"), and the amplitude of the largest negative peak (hereinafter 
defined more completely and referred to as "NVE"). 

The electrogram signal is filtered, preferably with a passbapd of two to thirty Hz (- 3dB). The lowpass sec- 
tion of the filter serves as an antialiasing filter, and smooths the transitions in the electrogram so that peaks 
may easily be identified. The highpass section of the filter serves to remove any baseline wander from the 

45 signal, so that a separate DC removal stage is not required. 

The filtered signal is sampled at a rate that can be as low as 80 Hz, using a known peak-holding method, 
such as the one discussed at pages 217-220 in the book The art of Electronics" (2nd Ed.), authored by P. 
Horowitz and W. Hill, published by Cambridge University Press, Cambridge, U.K. (1 989). Alternatively, the met- 
rics may be measured directly using analog electronics. The values of the metrics are calculated for a known 

so segment of normal sinus rhythm (NSR). The mean of the NSR value (hereinafter "NSR-mean") for each metric 
is stored for later use. 

When a rhythm is to be diagnosed, the NSR-mean is subtracted from each measured metric in the un- 
known rhythm. This gives the NSR-mean-removed metrics. If the unknown rhythm is in fact ST, then the values 
of all the NSR-mean removed metrics will be small, because the shape of the electrogram will be similar to 
55 that in NSR. If the rhythm is VT, however, then at least one (and usually more) of the NSR-mean-removed met- 
rics will be large. This difference is detected by the classification process, which is described below. 

The classification process includes a classifier, hereinafter described in greater detail and referred to as 
the "SUM classifier". The SUM classifier is an improvement over previously known methods of classification 
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theory. It provides similar performance at a much reduced computational load compared to previously known 
methods. In a step-wise fashion, the SUM classifier performs the following operations: 

a) Filters the electrogram; 

b) Calculates the raw metrics from the filtered electrogram; 

c) Amplitude-normalizes the metrics; 

d) Subtracts off the NSR-means; 

e) Subtracts off the ST-offeets; 

f) Takes the absolute values (to get error terms); 

g) Weights the error terms with weighing parameters; 

h) Sums the weighted error terms (to get the error sum); and 

i) Compares the error sum to a threshold. 

The foregoing steps give a classification of successive heartbeats as either normal or abnormal. A normal 
rhythm (i.e. p sinus tachycardia) is then defined as a rhythm that contains some number of normal heartbeats 
in a specified number of heartbeats, for example 7 or more normal heart-beats in the last 10 heartbeats. 

The metrics values are largely insensitive to filter corner frequencies. Sample rate also has little effect 
on performance, provided that peak-holding is employed. Good performance has been demonstrated at sam- 
ple rates as low as 80 Hz, as indicated earlier. 

Quantization of the signal to 8 bits causes only a small degradation in performance, and this is the rec- 
ommended signal resolution. It may be possible to quantize the signal further, but this has not been tested. 

The advantages of the method of the present invention for use in ICD's are, thus, that it employs a low 
sampling rate, requires very few multiplications per heartbeat consumes low power, and can be implemented 
in analog electronics. 

DESCRIPTION OF THE DRAWINGS 

While the specification concludes with claims particularly pointing out and distinctly claiming the subject 
matter of the invention, it is believed that the invention will be better understood from the following description, 
taken in conjunction with the accompanying drawings, In which: 

FIG. 1 shows the method of calculation of the MOD metric; 

FIG. 2 shows the method of calculation of the PKP metric; 

FIG. 3 shows the method of calculation of the DIF metric; and 

FIGS. 4A-4I show the values of the QRS metric for a number of representative wave shapes. 
DESCRIPTION OF THE PREFERRED EMBODIMENT 

A number of concepts recur in the following description 5 of the metrics. Referring to FIG. 1, which illus- 
trates a QRX complex 20 of an electrogram, the fiducial point or sensing point is shown at 10. The fiducial 
point 1 0 is a timing reference within the electrogram. In the simplest case the fiducial point is the point at which 
a sensing circuit (not shown) is tripped. The fiducial point that is currently used is the time of greatest dV/dt 
(slew) of the electrogram or, equivalent^, the time at which the high-pass-filtered electrogram has the greatest 
excursion from a baseline 12. 

An analysis window 14 is defined in terms of absolute offsets 16 and 18 from the fiducial point 10. In tests 
that have been performed utilizing the present method, the analysis window 14 has been set at 100 ms wide, 
centered on the fiducial point 10. The width of the analysis window 14 is not critical, but should extend out to 
the isoelectric segments on either side of the QRS complex 20 under consideration. Incoming (raw) data is 
filtered in an analysis filter (not shown) before any metric analysis is performed on the incoming data stream. 
The analysis filter removes low frequency base line wander, and high frequency noise. It also serves as an 
antialiasing filter. The analysis filter is preferably a first order 2 Hz highpass filter, followed by a second order 
30 Hz Butterworth low pass filter. 

In the discussions which follow, references to an electrogram are intended to mean the electrogram at the 
output of the analysis filter, with the exception of the metric hereinafter referred to as "FIL" metric, which metric 
uses a bandpass filter (not shown, but hereinafter discussed) in place of the 30 Hz Butterworth low pass filter. 

The following metrics which may be used in connection with the present invention will now be described 
in greater detail: MOD, PK, PVE, NVE, PKP, DIF, FIL, SSA, QRS and ORD. 

MOD Metric 



The MOD metric is the absolute (modulus) area under the QRS complex 20 of FIG. 1. It is shown shaded 
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at 22 and 24 in that figure. The MOD metric is defined as the SUM of absolute values of the samples within 
the analysis window 14. For an analogue implementation, the MOD metric would be defined in terms of an 
integral. This metric is intended to reflect the width of the complex 20 of the electrogram. 

5 PK Metric 

The PK metric is the value of the largest excursion of the electrogram complex 20 from the base line 12. 
It is signed so that if the largest peak is negative, than the PK metric is negative, and if the largest peak is 
positive (see, e.g. deflection 28 in FIG. 3), then the PK metric is positive. It reflects both the amplitude and 
10 polarity of the electrogram complex 20. 

PVE and NVE metrics 

The PVE metric is the amplitude of the largest positive deflection (see, e.g. deflection 28 in FIG. 3) in the 
15 electrogram complex 20. Similarly, the NVE metric is the amplitude of the largest negative deflection (see, 
e.g., deflection 30 in FIG. 3) in the electrogram complex 20. 

PKP metric 

20 The PKP metric is the peak-to-peak amplitude of the electrogram complex, and is shown in FIG. 2 at 26. 
This is the absolute amplitude of the largest positive deflection plus the absolute amplitude of the largest neg- 
ative deflection. It reflects the amplitude of the electrogram complex 20. 

DIF metric 

25 

Referring to FIG. 3, the DIF metric is the difference between the absolute amplitudes of the largest positive 
deflection 28 of the electrogram complex 20 and the largest negative deflection 30 of the electrogram complex 
20. The DIF metric will be large and positive for an upright electrogram complex, small for a biphasic one, and 
large and negative for an inverted one. Thus, the DIF metric is a continuous variable which reflects the polarity 
30 of the QRS complex of the electrogram. 

FIL metric 

The FIL metric is the peak-to-peak amplitude of a filtered version of the electrogram. The raw signal is 
35 passed through a filter having the following characteristics: single pole high pass at 2 Hz, two pole band pass, 
and Q = 1 at 25 Hz. The peak-to-peak amplitude of this signal is measured in the same way that the PKP metric 
is measured. This metric is designed to reflect the slew rate (dV/dt) of the electrogram complex 20. The high 
pass characteristic in the filter acts as a differentiator. 

40 SSA metric 

The SSA metric is the signed-square area metric. It is an integral measure, like the MOD metric. It is cal- 
culated as follows: for each sample within the analysis window the sample is squared and its sign is restored 
(e.g., 2 3 = 4; -2 3 = -4). Then the signed squares are summed, the square root of the absolute value of the sum 

45 is taken, and the sign is restored. This gives a measure having characteristics similar to those of the DIF metric, 
i.e., it forms a continuous variable that reflects the polarity of the electrogram. The squaring is used to ensure 
that a large amplitude narrow lobe in the electrogram will predominate over a smaller amplitude but wider one. 
This is necessary because the high pass characteristic in the analysis filter will often produce a wide, low am- 
plitude lobe in the filtered electrogram. The following is a mathematical description of the SSA metric: 

50 t = Z(sign(x).x 3 ) 

SSA=sign(t).VTtT 

QRS metric 

55 

Referring to FIGS 4A-4I, the QRS metric represents the QRS angle / as calculated from the amplitude and 
temporal order of the largest positive and largest negative peaks, as follows: 
r = abs(smaller peak larger peak) 
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H = arcos(r) 

If the first peak is positive and larger /=H 

If the first peak is positive and smaller: /=. H 

If the first peak is negative and larger /=H-1 80 

5 If the first peak is negative and smaller /=1 80-H 

ORD metric 

The ORD metric reflects the order of the peaks in the electrogram complex 20. It is used in conjunction 
10 with the DIF metric to replace the QRS metric. The ORD metric is defined as follows: 
If the positive peak is first, ORD = smaller peak + larger peak. 
If the negative peak is first ORD = -smaller peak + larger peak. 

This gives a metric (ORD) which is a continuous variable that reflects order. Thus, if the electrogram com- 
plex 20 is biphasic with the positive peak first, the result is close to one; if it is biphasic with the negative peak 
15 first, the result is close to minus one; if it is monophasic, the result is close to zero. 

The ORD and DIF metrics, in combination, would give ordered pairs of numbers that correspond very 
roughly to the QRS angles / of FIGS. 4A-4I, as shown in Table 1 below: 



Table 1 



- Relationship Among QRS, ORD 


and DIF Metrics 


QRS angle / 


(DIF.ORD) 


0 


(0.0,1.0) 


45 


(0.5,0.5) 


90 


(1.0,0.0) 


135 


(0.5,-0.5) 


180 (=-180) 


(0.0,-1.0) 


-135 


(0.5,-0.5) 


-90 


(-1.0,0.0) 


-45 


(-0.5,0.5) 



All three of the metrics QRS, ORD and DIF suffer from the problem that they represent a two-point de- 
40 scription of the electrogram complex 20, and some electrogram complexes necessarily require three points 
to properly describe them. They are simple to implement, however, requiring only positive and negative peak 
pickers (not shown), which are needed in any event for other metrics. 

As will be discussed in greater detail hereinafter, tests were performed to determine the correlation be- 
tween various parameters. All of the metrics, except QRS angle /, were highly correlated with each other. This 
45 is because they all depend very strongly on the amplitude of the electrogram complex 20. This does not rep- 
resent a problem for the PK and PKP metrics, which are intended to reflect the amplitude of the electrogram. 
However, it is unsatisfactory for the other metrics, which end up reflecting the amplitude more strongly than 
the desired feature. The result is a large collection of amplitude measures, and very little else. 

50 Amplitude Normalization 

To resolve the foregoing problem, metrics that have an unwanted sensitivity to amplitude are normalized 
in accordance with this invention. Such normalization is achieved by dividing the metrics in question by the 
PKP metric. This results in a set of metrics which are only loosely correlated with each other. Thus, they should 
55 reflect different aspects of the shape of the electrogram complex, as is desired. It is preferable to always use 
the MOD, PVE, NVE, DIF, FIL and SSA metrics in the normalized form. 

This process of amplitude normalization is preferably integrated into the SUM classifier, described below, 
enabling the number of computer operations to be minimized. 
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SUM classifier 



The SUM classif ierfor the metrics involves a complex process which decides if a particular measurement 
resides in a region of the measurement space that is considered normal (to be classed as ST). If the meas- 
5 urement falls outside of this region, it will classed as VT. Using the metrics MOD, PKP, and NVE f the meas- 
urement space will be three dimensional, and the normal region will be a three dimensional space around the 
mean of the ST measurements. 

The goal of the SUM classifier is to minimize the computational cost of its operation. This is primarily de- 
pendent on the number of multiplications made per heart-beat. 
10 Given a three dimensional feature vector, for example the three measurements MOD, PKP and NVE, de- 
noted m, p and n, respectively, and a set of weights (a, p, and y) the SUM classifier is defined as: 



ST 

_ _ < 

a\m-m ST \+V\p-p ST \+y\n-n ST \ T 

VT 



The above expression means that if the weighted sum of the distance of the measurement from the ST 
20 mean (ST offsets) is less than some threshold (T), the measurement will be classified as ST; otherwise, it would 
be classified as VT. Stated another way, it means that if the measurement is similar to the ST mean it will be 
classified as ST, and if it is dissimilar, it wilt be classified as VT. 

The foregoing discussion assumes that the measured metrics would be used directly. In fact, they must 
be processed somewhat before being used in the SUM classifier. Firstly, the m (MOD) and n (NVE) metrics 
25 need to be normalized by division by the p (PKP) metric, as follows: 

™' = f (2) 

30 The NSR-mean is then subtracted from the normalized metrics: 

mn=m '-™Ns* = J-m'NSR ( 4 ) 
P" = P-Pnsr (5) 

35 n-.n'-lk-jj-iL (6) 

where m ' 0TS , for example, is the mean of the normalized value of the MOD metric, calculated in NSR. 
The discriminant function is thus: 




45 This must be reorganized firstly to remove divisions, and then to minimize the number of multiplies that must 
be performed beat-by-beat: where 



50 



0 



m"nj, p"™, and n"™ are parameters (and 
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10 



15 



20 



25 



40 



ST 

^|*-p(ita + *jr> l*P-||p-tf«+P«J k^|ia-p(ii«*5«-) I < P < 8 > 



> 



so need not be calculated beat-by-beat), and 



^iTo Pi 



are calculated by background processes during NSR. This updating preferably occurs at a programmable In- 
terval In the range of 10 minutes to 1 hour. Thus, the number of multiplications required per complex during a 
rhythm classification is six. 

The terms a, 0 and y are calculated as the differences between the inverse of the square root of the di- 
30 agonal terms of the class covariance matrices K, described in greater detail below. 

The classifier parameters for use in an implant are calculated from a set of metrics produced by the real 
implant hardware and software. The parameters a, 0, y, m M TUf p'Vu. and n"™ are preferably programmable, 
but should not be available to the physician. This will allow the performance of the system to be tuned. The 
threshold parameter T (or an offset on the threshold) could be made available to the physician. This would 
35 allow a trade-off between sensitivity and specificity. 

In order to more clearly understand the term "covariance matrices K", referred to above, the following tu- 
torial material with respect to random variables is provided. 

The expectation or mean of a random vector V is a vector "m" defined as: 



w=E[x] =f~xp x (x)dx (12) 



where Py(x) is the probability density atx. This is taught by Therrien, C.W., in the book "Decision Esti- 
45 mation and Classification", published by John Wiley and Sons, New York, 1989, (hereinafter referred to as 
Therrien"), at equation 2.41 , page 24. 

When training a classifier from a limited data set, this is estimated as: 



50 




(13) 



where Xj k is the j'th observation vector for class 
55 i; Nj is the number of observations of class i; and 

mj is the estimate of the mean vector of class i. 
The covariance matrix "K" of a random vector "x" with expectation V is defined as: 

K=E[(x- m)(x-m)T] (14) 
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15 



20 



45 



where T is a matrix transposition operator (Therrien, page 49); and E is an expectation operator (Ther- 
rien, equation 12, page 54). (Also, see Therrien, equation 4.42, page 55). 
When training a classifier from a limited data set, this is estimated as: 

v- = JL ( ~ -A \ _A \T (15) 



10 where Kj is the estimate of the covariance matrix of class i. 

The SUM classifier relies on a specific discriminant function. Accordingly, the following discussion, relating 
to decision rules and discriminant functions, is presented at this time. 
The decision rule for for the two-class case takes the form: 



£{9) < T (16) 
> 



The interpretation of this expression is as follows. If some function f() of the observation vector y Is less 
than some threshold T, then assign the class label g>2 to it; otherwise assign the class label ©3. 

The SUM classifier is derived from a Gaussian quadratic classifier, which is well known in the literature. 
25 Accordingly, the following discussion relating to the Gaussian quadratic classifier is provided. 

The Gaussian quadratic classifier provides a Bayes optimal solution if the class conditional probability 
density functions are Gaussian. If the density functions are not Gaussian, then the probability of error is not 
minimised. The resulting classifier provides a decision surface which is best matched to the second moment 
statistics. (See, Therrien, pages 96,97). 
30 A general two-class quadratic classifier is defined as follows: 

h (y) -? r A?+b T ?+c < T ( 17 ) 

35 > 

(See Therrien, equation 6.6, page 96). 

This resulting classifier will be Gaussian quadratic if the classifier parameters are defined as follows: 
40 A = K? - K" 1 (18) 

b = 2(K?m 2 -K-;mi) (19) 
\k I 

c = m]K?mi - m\K?m 2 + Inj^j- (20) 

The derivation of the SUM classifier, starting with the equation forthe discriminant function of the Gaussian 
quadratic classifier, is set forth below. It develops the equations to calculate the parameters of the sum clas- 
sifier, viz. alpha, beta, and gamma. 

The discriminant function forthe SUM classifier can be derived from the Gaussian quadratic discriminant 
function, which can be written as: 

so h(y) = (y-mjrKffr-mi) - (y - m^K;^ - nij) +ln|^[ (21) 

In practice this will be compared to a threshold (for a two-class problem): 



55 
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h(y) T 
> 



(22) 



10 



15 



Thus the log of the covariance determinants can be considered as part of the threshold value T, resulting 
in the simpler form: 

h'{y) = (y - m A YK-;iy - #»,) - (y - m^ K ^(y - (23) 
If we set the class 2 (VT) mean to be equal to the class 1 (ST) mean, we can make the substitution: 

y- = y-m 1 (24) 

The discriminant function thus becomes: 

h"W =/ T (^ , -^V (25) 
If we set the off-diagonal terms for the covariance matrices to be zero, this becomes (for three features): 



20 






1 


0 


0 




1 


0 


0 




h /u {y f ) =y n 




0 


1 


0 




0 


1 


0 


25 






0 


0 


1 

o 2 cl 




0 


0 


1 

o\ 2 



(26) 



which reduces to: 



30 



35 



h'"(y') «jr* 



<?bl °b2 



(27) 



40 



which reduces to: 



45 



h l »(y') -y\ 



(28) 



50 



55 



From the above we see that the discriminant function is a weighted sum-of-square distance measure. (It 
actually measures the distances from the ST class mean, which intuitively makes sense.) 

In the special case where: < o£ < c&, and oi < this can be re-cast using a sum-of-absolute- 
values distance measure, giving: 

AGO ■ aly'J +p|y'J + y |y'J (29) 
where the coefficients a, 0, and y of the SUM classifier are thus calculated from the diagonal terms of the 
class covariance matrices as follows: 

a = V3TTjL (30) 



(31) 



10 



EP 0 653 224 A2 



el ^c2 

In practice it has been found sufficient to approximate this as follows: 

« = -7- (33) 

P = - L - - 1 - 04) 

T = ~ " 4- (35) 

It will be apparentf rom the foregoing description that the present invention provides a method of classifying 
heart rhythms by means of electrogram morphology, which method may be used in implantable cardioverter 
defibrillators, involves a low sample rate, very few multiplications per heart beat, can be implemented in digital 
as well as analog electronics, and consumes very little power. 

While particular embodiments of this invention have been shown and described, it will be obvious to those 
skilled in the art that various other changes and modifications may be made without departing from this in- 
vention in its broader aspects. For example, the MOD, PKP and NVE metrics are not the only suitable set of 
metrics that may be employed with this invention. Many other combinations of the metrics described herein 
may also be used. Similarly, other metrics may be devised that might also work well. In addition, the parameters 
of the SUM classifier could be calculated using equations 30-32, rather than 33-35. Alternatively, any other 
type of classifier known to the art could be substituted for the SUM classifier. These include discriminant func- 
tion, K-nearest neighbors, neural network, ...ID3. It is, therefore, aimed in the appended claims to cover all 
such changes and modifications as fall within the true spirit and scope of this invention. 



Claims 

1 . A method of classifying heart rhythms of a patient by means of electrogram morphology, comprising the 
steps of: 

storing a first set of at least two morphology-defining metrics derived from a QRS complex of a 
normal electrogram, said electrogram having a width and an amplitude, said metrics defining the shape 
of said electrogram including at least one of said width and amplitude ; 

establishing bounds of normal variation for said metrics which are combined to define a decision 
surface that includes first and second regions such that said first region represents normal heartbeats 
and said second region represents abnormal heartbeats; 

providing an electrogram of a patient's heart rhythm; 

measuring a second set of metrics for a given heartbeat from said patient's electrogram, said met- 
rics of said second set corresponding to said metrics of said first set; 

comparing the metrics of said second set with the metrics of said first set to decide which of said 
regions said given heartbeat falls into and to classify said given heartbeat as either normal or abnormal 
based on said decision; 

repeating said measuring and comparing steps for a plurality of times; and, 

if a number of heartbeats are found to be abnormal, diagnosing the patient* s heart rhythm as ven- 
tricular arrhythmia. 

2. A method according to claim 1 , wherein said storing step includes a sub-step of storing a first set of at 
least three metrics selected from a group of metrics that includes: a MOD metric, which comprises an 
absolute area under a QRS complex of a normal electrogram; a PK metric, which comprises a value of a 
largest excursion of said QRS complex from its baseline, said value having the same polarity as said ex- 
cursion; a PVE metric, which comprises an amplitude of a largest positive deflection of said QRS complex; 
a NVE metric, which comprises an amplitude of a largest negative deflection of said QRS complex; a PKP 
metric, which comprises a peak-to-peak amplitude of said QRS complex; a DIF metric, which comprises 
the difference between the absolute amplitudes of the largest positive deflection of said QRS complex 
and the largest negative deflection of said QRS complex; a FIL metric, which comprises a peak-to-peak 
amplitude of said QRS complex after said complex has passed through a bandpass filter having a center 
frequency of about 25 Hz; a SSA metric, which comprises a signed-square area metric in which a plurality 
of samples taken of said QRS complex are squared, their signs are restored, the signed squares are sum- 
med, the square root of the absolute value of the sum is taken, and the sign is restored thereto; a QRS 
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metric, which comprises a QRS angle 0 calculated from the amplitude and temporal order of the largest 
positive and largest negative peaks of the QRS complex, as follows: 

r = abs (smaller peak) + (larger peak), 

<D = arcos (r), and 

if the first peak is positive and larger, 0 = O, 
if the first peak is positive and smaller, 0 = -<D, 
if the first peak is negative and larger, 0 = <B - 180, and, 
if the first peak is negative and smaller, 0 = 180-0; 
and, an ORD metric, which is defined as follows: 

if the positive peak is first, ORD = smaller peak + larger peak, and 
if the negative peak is first, ORD = - smaller peak -!- larger peak. 

3. A method according to claim 2, wherein said storing step includes a sub-step of storing a first set of at 
least three metrics that includes said MOD metric, said PKP metric and said IMVE metric period. 

4. A method according to claim 2, wherein, if one or more of the MOD, PVE, NVE, DIF, FIL and SSA metrics 
are selected from said group of metrics, said storing step includes a further sub-step of amplitude nor- 
malizing said one or more metrics by dividing the value thereof by the value of said PKP metric. 

5. A method according to claim 1, including further steps, prior to said measuring and comparing steps, of: 

filtering the patient's electrogram through a bandpass filter having a center frequency of about 25 
Hz; and 

sampling said filtered electrogram at a predetermined rate to obtain said metric. 

6. A method according to any one of claims 1-5 including further steps, prior to said measuring and com- 
paring steps, of : 

said filtering of the patient's electrogram is performed through a filter having a pass band of about 
2to30Hz;and 

sampling said filtered electrogram at a predetermined rate . 

7. A method according to claim 6, wherein said sampling step includes storing a peak value of said electro- 
gram, and wherein said predetermined rate is about 80 Hz. 

8. A method according to claim 6, wherein said storing step includes a sub-step of: 

calculating and storing an NSR-mean value for each of the metrics of said first set, said NS R-mean 
value for each metric comprising the mean value of the metric during normal sinus rhythm of a normal 
electrogram. 

9. Amethod according to claim 8, wherein said comparing step includes a sub-step of subtracting from each 
of the metrics in said second set of metrics the corresponding NSR-mean value therefor to provide NSR- 
mean-removed metrics corresponding to the patient's heart rhythm. 

10. A method according to claim 9, wherein said comparing step includes a further sub-step of classifying a 
heartbeat as either normal or abnormal based on the value of one or more of said NSR-mean-removed 
metrics exceeding a predetermined threshold value. 

11. A method according to claim 10, wherein said diagnosing step includes a sub-step of diagnosing the pa- 
tient's heart rhythm as ventricular fibrillation if at least four or more heartbeats of ten successive heart- 
beats are classified as abnormal. 

12. A method according to any one of claims 1-11 , wherein said diagnosing step includes a sub-step of diag- 
nosing the patient's heart rhythm as sinus tachycardia if at least a specified number of a predetermined 
greater number of last heartbeats are classified as normal. 

13. A method according to any one of claims 1-11, wherein said diagnosing step includes a sub-step of diag- 
nosing the patient's heart rhythm as sinus tachycardia if at least least 7 of the last 1 0 heartbeats are clas- 
sified as normal. 
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